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Abstract. Wc prove via explicitly constructed initial data that solutions to the gravity- 
capillary wave system in MP representing a 2d air-water interface immediately fail to be 
with respect to the initial data if the initial data (/iOiV'o) € H^~^^ (g) i/* for s < 3. 
Similar results hold in domains with a Id interface. Furthermore, wc discuss the 
related threshold for the pure gravity water wave system. 



1. Introduction 



The motion of a perfect, incompressible and irrotational fluid under the influence of 
gravity and surface tension is described by the free surface Euler (or water-waves) equations. 
In this paper we consider the fluid domain to be given by fij = {{x , y , z) E Mj^ : z < 
h{t,x,y)}, where we have assumed that the the free surface is described by the graph 
z = h{t,x,y). The incompressibility and irrotationality of the flow indicate the existence 
of a velocity potential 0, and hence the Euler equation can be reduced to an equation on 
the free surface ([10], [11], [33]). More precisely, denoting the trace of the velocity potential 
(f) on the free surface by ip(t, x, y) = (pit, x, y, h{t, x, y)), the equations for ip and h are 



:i.i) 



dtip = —gh + T div 




2 1^*1 + 2(l + |V;f) ^ 
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with 

(Oxv)= (^') (xy)G^ ^^^^(^') ® ^^(^') ^hen r > 
where the (rescaled) Dirichlet- Neumann operator is given by 

and n+ is the outward unit normal vector to the free surface. 

There has been an extensive study of the well-posedness of the water wave problem. 
Early works go back to Nalimov [22], Yosihara [32] and Craig [9] where Id free surface 
are considered, surface tension is neglected and the motion of the free surface is a small 
perturbations of still water. The general local wellposedness of 2d full water wave problem 
was established by Wu [2H|, see also Ambrose-Masmoudi [2]. The large time well-posedness 
is established in [50] . 

In the case of 3d gravity water waves, Wu ^Hj proved it is local wellposedness. Lannes 
[20] considered the same problem for finite depth. The global wellposedness was proved by 
Germain-Masmoudi-Shatah [18] and Wu [31] independently. 

For gravity-capillary waves, Ambrose-Masmoudi [3] considered the zero surface tension 
limit and provided a proof of wellposedness in the 3d case. Alazard-Burq-Zuily [Tj studied 
the regularity of the local solutions using a paradifferential formulation in the 2d and 3d 
cases. See also the work of Christianson-Hur-Staffilani [2], where along with the work 
of Alazard-Burq-Zuily [T], smoothing and Strichartz estimates were established for the Id 
interface problem. 

When the fluid is rotational and without surface tension, the problem has been shown 
to be wellposed by Iguchi-Tanaka-Tani [19] in 2d and Lindblad [21] and Zhang-Zhang [31] 
in 3d. Including surface tension, the wellposedness is proved by Ogawa-Tani [25l [26] in 2d 
and Coutand-ShkoUer [15] and Shatah-Zeng [27] in 3d. 

The water wave system (11. II) generates a substantial and rich family of equations that 
describe solutions in various asymptotic regimes, including KdV, Camassa-Holm, BBM, 
among others. These model equations can accept remarkably low-regularity initial data, 
and so it is natural to ask whether the original system can also accept low-regularity data. 
In this paper, rather than continuing to investigate the wellposedness of water waves, we 
look at on the onset of breakdown of regularity of the solution operator, or flow map. It is 
shown that equations (11.11) start failing to behave in a "reasonable" way in rather regular 
spaces, counter to the behavior of many of the approximate equations listed above. In 
particular we identify the value of the regularity parameter s for which the model fails 
to be solvable by standard iterative methods. Our result joins up with the wellposedness 
results of Alazard-Burq-Zuily [T] , see Theorem [T] 

Unlike the semilinear dispersive system, where in many cases a contraction principle can 
be applied even in a very low-regularity regime to obtain wellposedness, and consequently 
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the smoothness of the solution map, the water wave problem is by nature quasilinear, 
suggesting that the initial profile can determine the existence time and the regularity of 
the solution map rather strongly. On the other hand, we will still use a perturbative 
approach by considering the nonlinear part of the system as a perturbation of the linear 
flow and show that the corresponding solution operator experiences a severe singularity in 
certain Sobolev spaces near the zero solution based upon those nonlinear interactions. 

We prove, via explicitly constructed initial data, that the solution map for the gravity- 
capillary wave system in (representing a 2c? air-water interface) immediately fails to be 
with respect to the initial data if the data {ho^ipo) is in H'+2 ® for s < 3. This 
is primarily due to the influence of surface tension. Similar results hold for a Id air-water 
interface in M?. The work, while not addressing fully the issue of ill-posedness of the 
gravity-capillary problem, follows from prior results related to the study of ill-posedness 
by Bona-Tzvetkov |5]. The idea stems from looking at the second and third order iterates 
in a Picard iteration about the chosen initial data, which leads to a high-to-high frequency 
interaction: a solution starts off initially with small energy and Fourier transform supported 
primarily at high frequencies, but quickly generates a large energy at high frequencies. In 
contrast to the high-to-low frequency cascade as in [H |T3] for the study of ill-posedness 
in Schrodinger equations, where the Sobolev index is negative, our high-high frequency 
interaction is chosen because of the ill-posedness into Sobolev spaces of positive indices. 
The initial data will be described in Section |3l but it is chosen to maximize to the greatest 
extent possible the frequency interactions in the nonlinear terms. A ill-posedness result 
was also obtained in the study of the nonlinear Schrodinger equation on the sphere in [TJ E] , 
in which the dispersive estimates are known to have a loss of derivatives due to the lack of 
dispersion. 

We note here that since we are looking at an immediate failure of the solution map to 
be with respect to the initial data, the effects here are related very much to nonlinear 
interactions through the Dirichlet-Neumann map and surface tension. The reason we go 
up to for a breakdown of regularity is due to the regularity restriction of the Dirichlet- 
Neumann map - h needs to be at least Lipschitz. But the second iteration only provides 
breakdown of regularity at the threshold at 

It is natural to ask if this this breakdown of regularity of the solution map is in fact a 
failure of uniform continuity of the solution map (a stronger measure of ill-posedness), which 
has been studied in semilinear dispersive problems, see for instance Bourgain-Pavlovic 
[B], Christ- CoUiander-Tao [12] and Bejenaru-Tao However, because of the quasilinear 
nature, it is very challenging to fit our problem in the framework of |1] , where they provide 
a general method for C'^ ill-posedness of semilinear problems. Implementation of this 
framework in (11. ip (and hence proving that our obstruction is indeed C°-ill-posed) may be 
achievable if the behavior of the Taylor expansion of the nonlinear terms in (II. ip (which 
we call N , see equation (I1.3P below) can be very well understood at all orders. 
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Proofs of local well-posedness for generic quasilinear equations typically involve parad- 
ifferential iteration schemes with dependence on derivatives of the solution in the inhomo- 
geneity of each iteration as well as the iterated linearization of the evolution map. Hence 
one should not expect the solution map to be smooth in relation to a particular Sobolev 
space W^, but rather in terms of scales of those spaces such that the inhomogeneous terms 
can be smoothly evolved. Indeed, given a solution m to a quasilinear equation, linearizing 
about such a solution may involve terms of the form d^u, which should live in a less regular 
Sobolev space. In addition, the water wave problem is implicitly strongly hyperbolic in 
nature, meaning that one should not necessarily expect sufficient gain of regularity in the 
evolution through dispersive methods to make a smooth flow map possible. Hence, for a 
full ill-posedness result, one would need to control the flow map at all relevant scales for 
the equation, which we do not do here. The result here distinctly uses the implicitly linear 
structure of the underlying evolution operator and only takes such quasilinear effects into 
consideration in the inhomogeneous, nonlinear iteration terms. 

We note that a key part of our analysis assumes the existence of a natural polynomial 
expansion of the Dirichlet-Neumann map G{h). This is by no means a trivial assumption 
and requires a fair bit of regularity to hold to a certain order as made clear by the increasing 
number of derivatives appearing in the expansion of G{h), see for instance |16]. Using the 
results of Nicholls-Reitich |23l |23] on the analyticity properties of the Dirichlet-Neumann 
map, one can be assured of a valid expansion up to order k in 2d if we require h G H'^ for 
s > 2 with 

\\Gk{hm\Hs<\m^.-.. 

Such requirements are strong, but our minimal regularity threshold moves well beyond this 
particular regularity requirement. 

Finally, we remark that when we consider the pure gravity water wave system (r = 0), 
the construction indicates the solution map is not for data (/icV'o) e H'-2 ® H' for 
s < |. However, in order to make sense of the Dirichlet-Neumann map, one needs to 
consider, at the minimum, an initial height ho in the class of Lipschitz functions, and the 
Hq G iJ* for s < 2 threshold fails to be Lipschitz. Iterating the pure gravity system further 
to third order does not move past the formal | critical Sobolev exponent threshold. Thus, 
this result can only be considered a formal failure of the iteration process. 

1.1. Formal argument for onset of regularity breakdown. Though the problem is 
fundamentally quasilinear, one can simply ask if the breakdown of flow-map regularity is 
dominated by the second-order terms in the nonlinear operator A/". Expanding the pure 
surface tension to second order results in a problem that can be looked at via scaling. In 
particular we have the model 

ht = \D\t/j - divihV^lj) - iDKhlDli/j) 

(1-2) 1 o 1 o 

i^, = Ah + -{\D\^pf--\ViJf 
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where we expanded the Dirichlet-to- Neumann operator to second order. If we look for scale 
invariant solution {h,4') ^ {hx,ipx) then 

2 2 

hx = X"^h{\t, \3x) 

1 2 



)H'^ with s < §• As we will see, 



A short calculation shows blowup occurs formally for H^'^'-: 
this argument is insufficient to understand the obstacle to local well-posedness in either the 
gravity-capillary or pure surface tension problem, since the mean-curvature operator does 
not show up until third order, and at that point the associated pseudodifferential operator 
becomes dominant. 

Remark 1.1. We note that a similar expansion of the pure gravity water wave equations 
leads to the model system 



ht 



\D\^ - div(Wz^) - iDKhlDlijj) 



2 



IZ^I 



with scale invariant solutions {h,ip) i— )■ {hx,ipx) satisfying 

hx = X~^h{Xt,X^x) 
^x = X-''ij{Xt, A^x), 



which blow up for data in H^~^ ® for s < |. Hence one could conjecture | to be 
the lower regularity threshold of ( II. ip . see Section |6l In this case the model system does 
accurately capture the essential dominant problem to low regularity solutions. We note, 
however, that Hq G H'^ for s < 2 fails to be Lipschitz. 

1.2. Abstract C'^ regularity breakdown of flow map argument. We follow the 
iterative method of Bona-Tzvetkov [5], used to study the BBM equation. Extracting the 
linear operator in (11.11) we write it as 



( dth= \D\iP + [G{h)iP 
dt^ = irA-g)h^ 




Ah 



Setting 





rA-g 



{G{h)i) + Vh-Vi)f 
2(1 + |V/i|^) 



\D 
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we write our nonlinear system of equations as 



:i.3) 



where 



dt 



= C 

/i(a;,0) 
V'(x,0) 



■MM), 



ho 

^0 



Then the Duhamel formula implies 



:i.4) 



ho 



The linear flow operator is determined by 



L2H + Li^ 



where 



Li(e,t) = cos(A(|e|)t) 



L2(e,t) = sin(A(|e|)t) 
Ls{^,t) = -smi\{\m 



^ e and A(r) G {Age, ^st, A^,}, where 

(1.5) Agc(r-) = \/gr + rr^, 

1 3 

(1.6) A^t(r) = r2r2, 

(1.7) Ag(r)=^5ri 

In the following we will choose (7, r G {0, 1} and hence three possible wave functions. Our 
results differ distinctly between the cases r = and r 7^ 0. 

Following the ideas of [S] and [1] we take 



h{x,0) 
^(x,0) 



a 



ho 



where a > and ho and ifjo are functions specified below. Of course, when a = the solution 
of fll.ip is the trivial solution. Our goal is to show that the map (a, /iq, V^o) — ^ {hit),ip{t)) 
is not with respect to a (when a = 0). This in turn implies the solution map to to fll.ip 
is not (see below). 
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In Duhamel form, we write (11. 3p as 



{a,t) = e'^a ( ) + / e^'-''^''MMdt' . 



Here have abused notation and written A/'(q!, t') in place of Af{h{a, t'),ip{(y, t')) in order to 
make explicit the dependence on a. 

From assumed wellposedness we have 



(o,t) 



0^ 

o; • 



We can continue computing higher variations. At first order we have 



;i.9) 



9„ 



{a,t) 



jc / ''-0 



0=0 



which follows from Duhamel's formula 



da 



'0/ Jo 



and the expansion in Section |2] below. For the second iterate we have 



;i.io) 



dl 



=0 Jo 



Finally, for the third iterate we have 



1.11) 



(a,t) 



a=0 JO 



If our solution map is indeed C with respect to the initial data and the topology X*(M ) 



then necessarily 
(1.12) 



a=0 



< c 



ho 

^0 



for a uniform constant C. Therefore, to show our mapping fails to be C^, we need to show 
that for any C, there exists explicit initial data Hq, ipQ such that 

k 



:i.i3) 



Qk 



Q = 



> C 



ho 



X" 



The following result identifies when the solution map fails continuity at the origin for 
gravity-capillary, pure surface tension, and pure gravity water wave equations with a one 
or two dimensional interface. 

Theorem 1. Let = H^'^^ and assume r > and g G [0, oo). For any T > the 
flow map associated to (11. ip is not of class 
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• from X' to C([0, T];X') fors<3 when d = 2. 

• from X' to C([0, T]-X') fors<\ when d=l. 

Remark 1.2. Alazard-Burq-Zuily [1] proved that (11. ip with r > is well-posed in X* for 
all s > 3. Therefore, our result provides in some sense an endpoint to local well-posedness. 

Remark 1.3. Formally, we find from our construction that the flow map associated to 
f lTT]) is not (or C^) from Y' = H'^l ® H' to C([0,T]; Y') for s < | when r = and 
^ > and d = 2. 

In order to prove this result we will expand the operator to third order. This produces 
cubic terms in the Duhamel term involving the linear propagators. Restricting our non- 
trivial data to ipQ with support on a high-frequency sector produces a class of data that 
fails continuity. 



2. Expansion of M and /c-iterates of (II. ip 

In this section, we assume the free suface h G if*(M^) with s > 2 so that the Dirichlet- 
Neumann map G{h) has an analytic expansion (see Nicholls-Reitich [231 123]). The following 
proposition details the quadratic and cubic approximation of (II. ip in 2d. 

Proposition 2.1. We have 

(2.1) 5«Ar(a,t)U= fnV 



(2.2) dlU{a,t)l^^ 



where 



-2[div(/iiVV^i) + |Z^|(/ii|Z^|^i)]^ 

(IDiv^O'-lv^ip 



(2.3) [^^\ = e- 



and so from (ll.lOp 

(2.4) (^^J=£e(*-*')^9^Ar(0,tOrft'. 

Finally, we can define the third iterate h-^, tp^ in terms of the first and second iterate. In 
particular we have 



(2.5) dlX{a,t)l__ 



'dlAr,{a,t)\^^, 
dlAr2{a,t)l^,^ 
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where 

diNM.t)L=^ = -3 [div(/iiV^2) + \D\ (/ii|Z^|^2)] 

(2.6) - 3 [div (/isV^Ai) + \D\ {h2\D\ilj{)] 

+ 3 [A {hl\D\i^,) + \D\ (/i^A^i) + 2\D\ {h,\D\ {h^\D\^P,))\ 

and 

^l^f2{a, t) = -3r div {Vh,\Wh,\^) 

(2.7) -6|D|7/>i(/iiAV'i + |D|(/ii|Z}|^i)) 

+ 3([|D|V'i][|Z}|^2]-V^i-VV^2). 

We prove this in the following argument. The Dirichlet-to-Neumann map is denoted 

G{h)ij = Goi' + Gi{h)ij + G2{h)i> + .... 

From [21], we have 

G{0) = Go = \Dl 

where D = —iVx and further from ([21], equation (21) on page 114) we have the recursion 
relation: 



Thus, 



G,m =D (hD^) - ^ J^GK/i)^) 

(2.8) '=0 ^ / 

(/iD^) - \D\ (hGoih)^) 

=D iJiDi)) - \D\ {h\D\^l)). 

Using the definition of D\ 

(2.9) Gi{h)^l) = - div (hVifj) - \D\ {h\D\ilj) . 
The recursion formula with n = 2 gives: 



(2.10) 



G,ih)ij=\D\D(^^D^'^-^\Dr^(^ 







(2-/)! 



Gi{h)iP 



\D\D -D^ - \D\' -Go{h)^ + \D\ (hGiih)^) . 
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Using the formulae for Gi and Gq then rearranging terms: 

= ^\D\D (h'D^) - QlZ^p {h'\D\^) + \D\ {h [- div (W^) - \D\ ih\D\m 

(2.11) =-l\D\' {h'\D\i;) + l\D\D [h^Di,) + \D\ {h [div (WV^) + \D\ {h\Dm} 
= - ^\D\^ {h^\D\i>) + \D\ S^^Dih^Dtl:) + hdiv {hV^j) + h\D\ {h\D\ilj)^ 

If we use D = — zV in the first term in the curly braces: 
(2.12) 

G2{h)ij = - ^\D\^ {h^\D\ij) + \D\ |-i div{h^Vij) + hdw {hViP) + h\D\ {h\D\iP) 

Then the product rule gives: 
-- divQi^Vi)) + h div {hVi)) = -- {2hVh ■ Vtp + h'^A'^j) + h'^Atp + hVh ■ Vip 

2 ^ 

And thus putting this into fl2.12p . we have 

G,ih)^ = -^\D\' {h'\D\^) + \D\ ^h'A^ + h\D\ ih\D\i;)^ . 
Using the relationship A = — |-Dp gives: 

G,ih)4^ = ^-A {h'm) + ^\D\ {h'Ai^) + \D\ ih\D\ ih\Dm . 

2.1. The a derivatives of A/". We expand the solution {h,ip) in powers of a: 

(2.13) h{a) = ahi + ^a^h2 + 0{a^) and V'(a) = a^i + ^a^?/;2 + O(a^). 

We need to compute the expansion of the nonlinearity Af in terms a. We want expressions 
for d^Af, n = l,2, 3. 

2.1.1. ForTVi- We have 

Ml := G{h)ij - Go^. 

Plugging in fl2.13p . we get: 

A^i = G{h{a))i){a) - Ga%l){a) 

(2.14) = Gi{ahi + a^h2/2 + 0{a^)){a%l)i + a^%l)2/2 + 0{a^)) 

+G2{ahi + a^h2/2 + 0(a2))(a^Ai + a^%lj2/2 + 0(a3)) + . . . . 
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Note that Gn{ah) = a"'Gn{(y)- Similarly for all n there is an n-multilinear operat 
Gn{hi, . . . , hn) such that 

Gn{h)=Gr,{K...,h). 

Using these two facts in fl2.14p and then collecting powers of a gives: 

MM) = a^G^ih)^^ + «3 f^Giih2)iJi + ^Giih)^2 + G2(/ii)^i] + 0{a^). 



So we have 



and 



daAfi = 0, 
dlM = -2 [div (hV^) + \D\ ih\D\^)] 



dlMi = 3G'i(/i2)^i + 3G'i(/ii)V^2 + 6G'2(/ii)^i. 
Or rather, if we use the above computations for Gi and G2'- 

dlUi = - 3 [div (/iiV^2) + \D\ {h^\D\ij2)] 
(2.15) - 3 [div {h^ViJi) + \D\ {h2\D\i,,)] 

+ 3 [A {h\\D\,Pr) + \D\ (/i?A^i) + 2\D\ {hr\D\ . 

2.1.2. For M- Now let us compute the expansion of M- Indeed, from fll.ip . we have 



T 



div 



Vh 



Ah 



1 + |V/i|2 
= I + 11 + III. 
Taylor expansion and substitution of (I2.13P gives: 

1 



{Gih)ij + Vh-Vijf 
2(l + |V/i|2) 



Using this and ( 12.13^ in / gives: 



1 - -\Vh\' + 0{\Vh\^) = 1 - a'-\Vh\' + 0{a^). 



(2.16) I = T 



div 



Vh 



1 + |V/iP 



T 



Ah 



div V/i 1 - a'-\Vh 



rdiv Vh -a^Vh 



Ah 



0{a' 
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The similar expansion for III is pretty easy: 

11 = -^|V^p = |v^^|2 _ «3lv^i ■ V^2 + 

To deal with III notice that ^HM shows that {G{h)^ + V/i ■ Vz/;)^ = 0{a^). Also, 
Taylor expansion gives: 

1 



l + 0{\VhY) = l + 0{a' 



So we have: 

Then, since G2{h)ip = 0{a^), we have, with f l2.13p 

(2.17) =^ (aGoA + ^a^Go^2 + + a^Vh ■ V^i^ + 0{a^) 



2 

If we insert the expressions for Gi and Gq we get 
(2.18) 

/// =a'i(|D|^i)2 + a%\D\i;,) Q|D|^2 - div(/iiVV^i) - \D\{h\D\i;,) + ■ 
+ 0(a^) 

Differentiation with respect to a and organizing terms gives: 

9„Ar2 = 0. 

Again we differentiate: 

^lA^2 = {\D\^p^r-\V^P,\^ 

giving (1221). 

The third derivative is: 

dlAr2 = -3rdiv(V/ii|V/iip) -3V^Ai ■ V^/'2 

+6(|D|^i) Q|D|^2 - /iiA^i - \D\ih\D\i;^)\ . 
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A quick reorganization of terms gives: 



3rdiv (V/ii|V/iip) 

-6|D|^i (/iiA^i + |D|(/ii|D|^i)) 

+ 3([|D|^i][|D|^2]-VV^i-V^A2), 



(2.19) 



which completes (12.51) . 

For Id interfaces, the expansion of daM {p{a))\a=o produces six terms instead of the 
twelve terms in the 2d case. 



3. The Quadratic Threshold for Regularity of the Flow Map 

We begin by studying the second-order iteration of the solution map. From fll.lOp and 
(12. ip . we have the following Duhamel terms 



(3.1) 




Let 



div(W^) - 1^1 {h\D\i)), 



i(iDi^)^-iiv^r 



Then, 




where 



m2{^,v) = ^{i^-v)-v + 1"^ -^ll^l) 



21^-^11^1 (cos^2 + 1) 



and 



(3.2) 



angle between ^ and 77, 
angle between — 77 and rj. 
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Expanding Q in terms of the linear propagators yields 



(3.3) 



,it-t')C 



for the implicitly defined Qj. Taking the Fourier Transform of the above second order 
iterations yields the following set of terms: 

Qi = t - t')Qi(e, t') + M^, t - t')Q^{i, t') 

t - v) - V, t'm -v) + - v, - v) 

+ j L2(e, t - t')m2i^, ri) (L3(e - V, t'M -v) + Li{^- V, t'm - v)) 

and 

X fL3(r/,t')/J(^7) + ^i(^,t'Mr/))c?r/. 



Of the 16 terms in the two Duhamel operators Qi and Q2, we look for a dominant lower 
bound term that violates an assumed upper bound. This lower bound will arise from the 
following construction: 



3.1. Candidate for Breakdown of Flow Map Regularity. For large and 6 small 
we set 



(3.4) 



/i = 0, 



N<\^\<2N^e\<5, 
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where ^ = |^|e*^; hence, the support of ip hes on a small arc far from the origin. Then our 
expression for Qi and Q2 reduces to the following four terms: 



(3.5) 





~t'\ 










-t'\ 


)m2(e,r/)£i(e 








-t'\ 


)mi(e,r/)L2(e 


-r/,t')2i(r/,t'Me 






-t'] 


)^2(e,^)2i(e 




— r])il){r])dri. 



We will show that D dominates the other terms and gives the required lower bound. 
Since mi and m2 have explicit dependence on the angle between ^ and 77, we first prove 
an auxiliary lemma that explains how the support of our constructed if) affects the support 
where — rfjipirj) is nontrivial. 

In the following we write A < B ii there exists a universal constant C such that A < CB. 
Likewise we write A > 5 if there exists a C with A > CB. Finally, write A B ii A < B 
and A>B. 

Lemma 3.1. If N is large and 6 is small then ip^^ — T])ilj{r]) has support in a set with 

(3.6) \^\ N and \e\ < 5 
where ^ = \C,\e^^ ■ Furthermore, we have 

(3.7) \m,i^,v)\<N'6', 

(3.8) |m2(e,r/)|~7V2_ 

Proof. Supose ^ — r] and r] both lie in the same sector. Then by the parallelogram law, ^ 
lies in that sector too, and |^| < 6. Thus, we also have |^i| < S and |^2| < 25. Therefore, 

I cos^i - 1| < 572, 
1 < I COS^2 + 1| < 2 

for 5 < f . 

Obviously we have |^| < 4:N. On the other hand, by the law of cosines |'CP = l'C~'7p + 
|?7p — 2|,^ — ri\\ri\ cos(7r — 62), so 

|^|2 > 2iV2 _ 2]v2 cos(7r - 25) = 2A^=^(1 + cos 25) > 2N^. 

Hence ~ iV and therefore, 

l"^i(e,^)l = lell^l|cos^i-l| <iW 
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and 
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l"^2(e,^)i = ^ie-^ll^llcos^2 + i|~A^'. 



□ 



(3.9) 



In the gravity-capillary case we have \gc{r) = y/r^ + r so when r 1, 



T 2 . 



Next, in order for cos(Agc('")^') > 2 ^oi r 1 and all < t' < t, we require that t <^ r 2; 
therefore, we choose 



(3.10) 



11, 3 

< t < with a > -. 

- - 100 iV« - 2 



The following proposition follows from Lemma [3.11 and fl3.9p -f l3TT0|) . 



Proposition 3.2. On the support of^jj{^ — r])^jj{r]) and for < t' < t < with a > 3/2 

we have the following hounds: 

I L3(e, t - t')m,L2{^ - V, t')L^{7], t') + Li(e, t - t')m2L,{^ - r/, t')L,{7], t') 

+L2(e, t - t')m2{^, r])Li{^ - r], t')L,{r], t') + Li(e, t - t')m,{^, r^)!^^^ - r/, t')2i(^, t') 

Proof. From (13.101) and Lemma 13.11 we know that for large 

m{t - t'l mt\ m - v)t', kvY < ^^ni-'^ < 



Therefore 



and hence 



l>L,{i,t-t')Xi{i-V.t')Xi{v.t') > \ 



lb 



On the other hand. 



^37721^2-^1 



<\sin{Xm-t'))\^m\ 17 



(3.11) 



6' 



\^-riMX'm-t')t' 



5' 



sin(A(e-r/)tO 



A(e - V) 



sin(A(0(t-t')) 



xm-t') 



sin (A(e - v)t') 



4ie-.iHA^(0(t-t')t'4(4 
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Similarly, 



L2m2LiLi 



<|sin(A(0(t-tO)l3||yl^-^ll^l 



\i - V\m\it - t') 



sin(A(0(t-t')) 



<\i-ri\m\it-t')<N'~\ 



LimiL2Li 



6^ 

<j\m 



sin(A(^-r/)t') 



5' 



sin(A(^-r/)t') 



and we obtain the claimed lower bound. 



□ 



3.2. Breakdown of Flow Map Regularity. We now give a result relating to a 
regularity threshold. Assume (/i, ip) are wellposed in X* = if^+a ^H'^ for a range of s G M. 
Then 



sup \\Q\\xs < 

0<t<T 



ho 



X" 



Taking {ho,ipo) as in (13. 4p . then 
(3.12) 



L2 



^1/2 



since the support lies on a set of area {N'^5). So 



(3.13) 



sup ||Q||x= < S. 

0<t<l 



Recall that 0<t< Y^-^fora> |. Note that from Lemma [XT] the region where Q is 
nontrivial lies in the high frequency sector of ^. In particular, consider the high frequency 
sector E = {2N < \^\ < AN, \9\ < S/2}, then we know \E\ ~ From Proposition O 
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we have 
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sup II Q\\xs > 2 sup 
o<t<r o<t<T 

> 2 sup 

0<i<T 



(0 



- ' Qidt' 



(0^ / * Q2dt' 



+ 2 sup 

L2(£;) 0<t<T 



(0^ / Q2rft' 









L2(i?) 



L2(i?) 

1 

r lOOAT" /" „ „ „ 

_ jY3-a-s^3/2^ 

where / = (r; = |r7|e*^ : - 10 < |r/| < + 10, |6'| < 5]. Taking 5 = iV"^' for some e > 
then 

sup IIQII^, > A^°+(5, 

0<t<T 



which violates (13.131) when 0<3 — a — e — s<3 — 3/2 — e — s = 3/2 — 3e — sor 

s < 3/2. 

As mentioned before, this threshold corresponds to solutions with h G H'^ for s < 2. 
Since h is not necessarily Lipschitz, the Dirichlet-Neumann operator may not make sense. 
We next iterate to the third level, where the surface tension operator appears. 

4. Cubic Terms 

We take our candidate for regularity breakdown and see how it behaves under the next 
iteration. For A^ large and 6 small we again set 



(4.1) 



h = 0, 
^= A^-(^+i)l 



N<\i\<2Ni-\9\<5, 



where ^ = |^|e*^; hence, the support of ijj lies on a small arc far from the origin. We set 

I = {?]= \7]\e'^ : A^ < |?7| < 2A^, |^| < 6}. 



From (11. lip . (12.61) and (12.71) . the third iterate now can be decomposed into a quadratic 
interaction between the first and second iterations and also a cubic interaction of the first 
iteration. To be more precise, we have 



(4.2) C 



{t')dt' 



= 2,Q + C, 
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where 

Qi = - div (/iiV^a) - \D\ {hi\D\tP2) - div (/laV^/^i) - \D\ (/lall^l^i) , 

(4.3) Ci =3 [A {hl\D\i>^) + \D\ {hlAiJ^) + 2|D|(/ii|D|(/ii|D|V;i))] , 
C2= -6[hiAij, + \D\{h\D\^i)]\D\ij, 

First, we prove bounds on the cubic interactions resulting from the cubic terms in the 
nonhnear expansion. Let pi correspond to the appropriate cubic pseudodifferential multi- 
pher leading to the cubic terms in Ci. Let p2,i, P2,2 correspond to the appropriate cubic 
pseudodifferential multipliers leading to the cubic terms (/ii, /ii, ipi), respectively (/ii, /ii, hi) 
in C2. Expanding, we have 

Pi{i,v,^) = -3|e||r/|(|e| + |r/|-2|e-H), 
P2,i{^^J^,v) = 6|r/| [|z/p - |^-r/||z/|] , 

P2,2(^, i^, v) = -9r [(^1 -Vi- ^i)^ir{l + (6 - ^72 - z^2)z^2'72] 
[(6 -rii- vi)vir]l + (^2 - ?72 - T^2y2vV\ 
-12r [(^1 - r^i - i^i)z/2?7i?72] • 

Lemma 4.1. // ^s large and 5 is small then tpi^S, — i] ~ ^)i^{.'n)i^W) support in a set 
with 

(4.4) 1^1 ~ N and \e\ < 5 
where = \C,\e^^ . 

Proof. From (14. ip we know that C.—rj—i', rj and u all lie in the same sector. By parallelogram 
law, 7] + u also lies in the sector. Because ^ = — r] — i^) + [r] + i^), then by parallelogram 
law, ^ lies in that sector too, and \9\ < 6. As for an estimate of |^|, obviously we have 
\^\<\^-V-'^\ + \v\ + W\ < 6iV- 

By law of cosines, we can find angles ai and 02 between ^ — r] — u and t] + u, rj and u 
respectively such that 

..rj. lOi lO I II I 

1^1 = Is ~ ''7 ~ '^l + 1^7 + i^l — 2|^ — ?7 — z/||?7 + i/| costti, 
\v + 1^1^ = 1^7^ + Wl"^ — 2|?7||i^| cosa2, 

so 

\r] + z/p > iV^ + iV^ - 2N^ cos(7r - 25) = 2N^ (1 + cos(25)) . 

Therefore, 

|^|2 > + 2N^ (1 + cos(25)) + 2ArV2(l + cos(25))cos(25) > 4N^ 
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and thus 
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2N < \i\ < QN. 



□ 



The following lemma provides a useful bound on the pseudodifferential operators arising 
in the cubic expansion. 

Lemma 4.2. If N is large and 5 is small then we have 

(4.5) \P2,i{^,v)\<N\ 

(4.6) \p2,2{^,v)\-N\ 

Proof. The bound on p2,i is immediate. Since — u = — t] — u) + t], we know that C, — 
also lies in the sector |^| < S and |^ — ~ A^. Then the estimate of pi can be easily 
obtained from Lemma 14. 1[ 



To finish the proof, denote 

^ — r] — u = rie*^^ , t] = r2e*''^ , 
then \(3j\ < S, N <rj < 2N, j = 1, 2, 3. Hence 

\P2,2\ ~ 9rir3rl cos (3i cos (3^ cos^ /32 + 9rir3rl sin (3i sin 3(3^ sin^ (32 + ^rir^rl cos (3i cos (3^ sin^ /32 
+ Srir^rl sin /3i sin /^s cos^ + ^^rir-^rl cos /3i cos /^a cos (32 sin /32 

□ 

Recall that 0<t<Y^-^fora>|. Note that from Lemma [4.11 the region where C is 
nontrivial lies in the high frequency sector of ^. In particular, consider the high frequency 
sector E = {2N < \^\ < 6N, \e\ < S}, then we know \E\ ~ A^^^. 

We are left with understanding the Duhamel term 



(4.7) 
Then, 



C 



Ci(t,0 

C2{t,0 



\{t-t')c ( Ci(t',/ii,V^i)\ , 

C2{t'M.^l) 



P2,i (^, J^, v)Ht, i-r]-v)^i (t, z/)^i (t, T]) 
+P2,2(C, ^:Tl)hi{t,i- T] - u)hi{t, u)hi{t,'r]) drjdu. 
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Expanding C in terms of the linear propagators yields 



21 




dt' 



for the implicitly defined Cj. Taking the Fourier Transform of the above second order 
iterations yield the following set of terms: 



implying 



X (L2{v,t')h{r]) + Li{r],t')ip{7])^ dr^dv 



and a similar expression for C2 



C2 = L2{i, t - t')C^{i, t') + Li(e, t - t')C2{i. t'), 
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which imphes 

4 = 




X [L^{i -r,- u,t')h{^ -r]-u) + Lgl^ - r] - u,t')^{^ - r] - u) 



+ J J Li{^,t- t')P2,i(e, ly, V) [H^, t')h{u) + Li(z/, t')V^(i^)^ 
X - r/ - I., t')h{^ -r]-u) + L;{^ - r/ - i/, - V - ^) 

As before, we may only look at €2- Setting = and collecting like terms, we obtain 



+ t - t')p2,i(e, r/)23(e -V-'^, t')Li{u, t')Li{r], t') 

t - t')P2,2(e, ^)23(e -n-u, t')L-siu, t')L^{r,, t')) ^^(e - ^ - y)^{u)^{n)dr,dv. 
Choosing < t' < t < Yoo^Wa' '^^ '^^^ bound 



Moreover, we have 

L2P1L3L3L1 + LiP2^iLiLiLi 
Since 



L2m,{u) = sin (A(0(t - t')) tItt ■ sin (A(i^)t') « 1, 



and by Lemma HI 

L2P1L3L3L1 + Lip2,iL3LiLi 



\i-ri-v\ 



sin (A(^ - - . 
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The third multipher in C2 can be estimated in the following way: 
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Li P2, 2-^3-^3-^ 



sin (A({ -T]- p)t') sin (A(z/)t') sin (A(r/)t' 



A(e-^-z/)A(z/)A(r/) 



1^ -?7 - v\\v\\ri\ 

sin (A(^ -r]-u)- X{v) + A(r/)) t' - sin (A(^ - r/ - z/) + A(z/) + A(r/)) t' 



iV"/2 gin (^(^ _ ^ _ j,)^/) gjj^ {\{v)t') sin (A(r/)t') 

iVll/2 



+ sin (A(^ -ri-u) + X{u) - A(r/)) t' - sin (A(^ - rf - u) - X{u) - A(r/)) t'] 



Therefore, 



sup 

0<t<T 



> 6 sup 

0<t<T 



> 6 sup 

0<t<T 

^ sup 

0<t<T 



^0 

Jo 

■t 



+ 6 sup 

L^{E) 0<t<T 







L2(B) 



(0^ 




-^1^*2,2-^3-^3-^3 

dridpdt' 

LHE) 



First, we have 



sup 

0<t<T 



(0^ 




L2P1L3L3LI + LIP2AL3LILI 



7572 



■j/''?/''?/^ dr]dvdt' 



L2(i?) 



L2(i?) 



7V1/2 



-2s 



i — COS 



'/J/ A(^-r/-z/) 

^l/2-2sjy-3/2|j|2|^|l/2 _ ^l/2-2s^-3/2^4^2^^1/2 
^4-2s^5/2_ 



L\E) 
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Similarly, 



sup 

0<t<T 




-^1^2,2-^^3-^3-^3 



d-qdudt' 



■2s 



1 - COS [(A(^ -7]-u)- X{u) + X{r])) Ar-3/2/ioo] 



_ ^5/2-2. jy-3/2|j|2|^|l/2 
_ ^6-2.^5/2^ 

where we have used the fact that in the support of ^'{^ ~ V ~ ^)i^{v)i^{^)y 

m --n-y)- m + m\ , m - ^ - ^) - m - \{ri)\ , 

|A(^ - 7/ - J.) + X{y) + A(r/)| , \X{i - - u) + A(z/) - A(r/)| ~ N'^'. 

Thus, 



drjdv 



(4.8) 



sup 

0<t<T 



C 



Now, we seek to prove bounds on the quadratic terms in the expression for d^f/. In 
particular, we have 



Qi 



fni{hitp2 + h2'4'i) dr], 



where Qi ^"^^ Q2 are given in (14.31) . 

As a result, the contribution to the cubic iteration from the quadratic terms is 

Q2 = L2Q1 + L1Q2, 



where 



2m ^ At-t')c ( Qi 
Q2J \Q2 
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Let us take Q2 and analyze its size in H^: 



Q2 



t PS 



L2mi{Liho + Ls'ilJo)Q2{t')dr)dt'ds 



Jo 

t PS 



^0 

t PS 



JO 



L2mi{L2ho + Li^lJo)Qi{t')dr]dt'ds 
Lim2{L2ho + LiilJo)Q2{t')dr]dtds, 



where we have used (12.21) and (12 ■4p . and Qj is defined in (13.31) . Taking < t < jq^p^k for 
a > 3/2, bounding Q2 in gives 



t PS 



iO" III L2mi ( L3i'oQ2 + LitpoQi ) + 2Lim2LiV^oQ2 cir/citcis 
Jo 



LHE) 



where we have used that 

L2miLs + L2miLi + Lim2Li 



< j N^N-'' dri ~ A^^A^-^'|J|. 



Therefore choosing a = 3/2 we see that 



Q2 
S 



< jY6-2a-2s^5/2 _ p^3-2s^5/2 



A similar result holds for 

For 2d interface problem we can now finish the 
Proof of Theorem Ui We know that 



Hence it is of lower order to the largest cubic term. 



(4.9) 



ho 

^0 



oWh- 



L2 



6i 



since the support lies on a set of area So 



(4.10) 



\C\\^<5'/'. 



sup rZ 

0<t<l 



Taking 6 = N from (14. 8 p we have 



sup ll-ll;^. ^ 

0<t<T 



which violates (I4.10p when 



6-2s-2e>0 
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or 

s < 3 - e. 

□ 

5. Id INTERFACE 

We make a few comments on the proof of the regularity threshold in Id. As in the 2d 
calculation, we restrict tpo to an interval at high frequency and set 

h = 

(5-1) 

As in the 2d case, the quadratic iteration provides a regularity threshold in = H^~^^ i^H^ 
for s < 3/2, which is borderline Lipschitz for h in one dimension as required to make sense 
of the expansion in [23] . Iterating further we find, as in the 2d calculation that 

(5.2) C= f e(*^*')^9;^Ar(0, t')dt' = 6 /* e^*"*')^ ( {t')dt' + f e^*"*')^ ( {t')dt' 



Q2J ' ' Jo \C' 



where 



(5.3) 



Qi=-d^ {hid^tp2) - \dx\ {hi\d:,\ip2) - (/i24^/'i) - \dx\ (/i2|9x|V^i) , 
Q2 = \dx\i^i\dx\ip2 - dxtpidxtp2, 

Ci =3 [dxx {hl\dx\4Ji) + \dx\ {hldxx^Pi) + 2\dx\{h\dx\{h\dx\tPi))] , 
C2 = - 6 [hdxxi^i + \dx\{h\dx\^i)] \dx\^i - 9r [dxxhidxhy] . 

As in the 2d case, the terms arising from the second-order iteration, Qj, are lower order, 
and we can restrict our attention to the terms arising as cubic interactions of the linear 
propagators hi,ijji. The cubic order pseudodifferential operator multipliers are simpler in 
the Id case, namely we have 

Pi = -m\v\m + \v\-m-i^\), 

P2,2 = -9r(^ -7]- v)v7f 



and ~ V ^ '^)^{v)^i'^) has support on the set E = {3N < ^ < 6A^}. As in Lemma [4.21 

Pi < N', 

P2,l < 
P2,2 ~ N"^- 
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We again concentrate on C2 which produces the largest term of the cubic expansion. Again 
we have 

+ Li(e, t - t')p2,i{^, u, r])Ls{^ -V-^, t')L^{u, t')Li{r], t') 

+Li(e, t - t')P2,2(e, r,)L^{i -V-^, t')L-siu, t')L^{r,, t')) - r] - u)ip{u)ip{r])dr]du. 
We study the third muhipher in C2, 

-^1^2,2-^3-^3-^3 , 

which is again bounded by 

iVll/2 



sin (A(e - r/ - z/) - A(z/) + A(r/)) t' - sin {X{^ - t] - u) - X{u) - A(r/)) t' 
'I 

- sin (A(e - 7/ - z/) + A(z/) + A(r/)) t' + sin (A(e - r/ - z/) + A(z/) - A(r/)) t'] . 
Thus, if / = {A^ < ^ < 2iV}, then 



sup 

0<t<T 



(0^ 



'0 

^4-2 




-^1^2,2-^3-^3-^3 



drjdvdt' 



1 - COS [(A(e - - z/) - A(z/) + A(r/)) Ar-3/2/iQo] 



iJi A(e-r/-z/)-A(z.) + A(r/) 

^4-2sjY-3/2|j|2|^|l/2 
5-2s 



drjdv 



L\E) 



where we have used the fact that in the support of '(/'(■C ~ V ~ ^)i'{v)'4'{^)j 

|A(e + m + A(r/)| , |A(e - - z.) + A(z/) - A(r/)| ~ A^^/^ 

As a result, 



(5.4) 

We know that 
(5.5) 



sup 

0<t<T 



C 



ho 



L2 



since the support lies on an interval of length A^. So 



(5.6) 



sup < 1. 

0<t<l 



28 R. M. CHEN, J.L. MARZUOLA, D. SPIRN, AND J. D. WRIGHT 

Then from (15.41) we have 

sup ||C||^. > » 1, 

0<t<T 

which violates (15. 6p when 5 — 2s > or 

s < 5/2. 



6. Pure gravity problem 



As mentioned in the introduction we visit here the flow map regularity threshold for the 
pure gravity problem. In this case, we gain nothing by going past the regularity threshold 
at the level of the quadratic terms, as the cubic interactions due to surface tension are 
responsible for the improved results above. As a result, we proceed as in Section [31 In the 
pure gravity case, we recall from (II. 7p that we have Xg{r) = ra. Hence, when r ^ 1, 



(6.1) 

and we choose 
(6.2) 



A.(r) 



r 2 . 



11, 1 

< t < with a > -. 

- - 100 iV^ - 2 



in order for cos{Xg{r)t') > | for r ^ 1 and all < t' < t. Following as above we can show 
that 

Proposition 6.1. Given Li, L2, L3 defined using Xg in (II. 7p . on the support ofil){^—r])il){r]) 
and for < t' < t < with a> 1/2 we have the following bounds: 

L^i^, t - t')m,L2{^ - f], t')L^{r], t') + L^{^, t - t')m^L^{i - r,, t')L,{r], t') ~ 



Proof. Note, this statement is exactly analogous to that of Proposition 13.21 The proof will 
follow quite similarly, except that we have now 

mit - 1% A(ot', A(e - v)t', mt' < ^ATi- < 1. 



Therefore, 



and hence 



L^{i,t-t')m2L^{i-r],t')L^{r],t') > —\^ - r]\\r]\ 

16 
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On the other hand 
L3miL2Li 



A(0 



<2|sin(A(0(t-t'))l^lell^lh7 



(6.3) 



s'\^-v\m\mt-ty 



(5^ 



sin(A(^ - r])t') 



sin(A(0(t-t')) 



<6'\^-vMym-t')t'<6 



1 



sin (A(^ - r])t') 



50 



So we get the claimed lower bound. 



□ 



Remark 6.1. Notice that in the above proposition we do not include the other two multi- 
pliers L2(e, t - t')m2i^, r/)Li(e - V, t')Li{ri, f) and L,{^, t - t')m,{^, 7])L2{i - r/, t')L^{r], t') 
like in Proposition 13.21 This is because in obtaining a lower bound of supQ<j<j. || Q||ya, 

we only need to estimate || supQ<j<7^ Q2dt'\\H'>-, which involves only two multipliers given 

in Proposition 16.11 The estimates for the other two multipliers L2m2LiLi and Lim\L2Li 
are used at the third iterates. However in the gravity case, to obtain similar estimates for 
-^27722-^1^1 and -Limi-L2-Li, one would need to require a > 1, which lowers the regularity 
threshold s. Therefore going up to the third iterates for the gravity waves does not help to 
improve the threshold exponent. 



For the pure gravity case we follow the same argument with a few changes. Assume 
(/i, ip) are wellposed in = H'^'^ (g) for a range of s G M. Then 

2 



sup 

0<f<T 



< 



ho 



Assuming initial data as in f l3.4p and the resulting bound fl3.12p . then 

sup IIQIIy. < 5. 

0<t<l 

Recall that < t < for a > 1/2. We again consider the high frequency sector E as 

before. From Proposition 16. II we have again 



sup 

0<t<T 



\Q\\yb > 2 sup 



0<t<T 



Qidt' 



L\E) 



2 sup 

0<i<T 



{iY / Q2dt' 



L\E) 



therefore, 



sup 

0<t<T 



violates f l3.12p so long as 

0<3-a-e-s<3-l/2-e-s 



5/2 



e — s 
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